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Real Corrections to the Virtual Photon Impact Factor 

J. Barteis*, S. Gieseke* 1 ', A. Kyrieleis 1 " 

II. Institut fur Theoretische Physik, Universitdt Hamburg 

Luruper Chaussee 14-9, 22761 Hamburg, Germany 
O ' 
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CN ' We calculate, for the longitudinally polarized virtual photon, the cross 

section of the process 7* + q — ► (qqg) + q at high energies with a large rapidity 
gap between the fragmentation system qqg and the other quark. This process 

C^ I provides the real corrections of the virtual photon impact factor in the next- 

to leading order. Evidence is given for the appearance of a new qqg Fock- 
component of the photon state. 



I. INTRODUCTION 



After the completion of the NLO virtual corrections |]J, the next step in our program of 
computing the virtual photon impact factor in next-to-leading order is the calculation of 
the emission of a real gluon, as sketched in Fig. [1|. We will compute these corrections from 
the high energy limit of the process 7* + q — > (qqg) + q (Fig. |2|) at leading order in a s . 
We are interested in the region of phase space, were the gg-pair and the g are emitted in 
the fragmentation region of the virtual photon with a large rapidity gap relative to the 
outgoing lower quark. We impose no further restrictions on the mass of the (ggg)-state, and 
therefore our results will also include the production of the gluon in the central region: this 
corresponds to the leading-order BFKL calculation, where the gluon is emitted with a large 
rapidity gap relative to both the lower quark and to the upper gg-pair. For the calculation 
of the photon impact factor, this piece has to be subtracted from our results. 

A study of the process 7* + q — > (qqg) + q includes the subprocess 7* + g — > qqg with 
an off-shell incoming gluon. Together with the NLO corrections to the production process 
7* + g — > qq which has been calculated in |I[, it provides the complete NLO corrections to 
the photon-gluon fusion process with off-shell incoming bosons. This generalizes the usual 
on-shell calculation of the collinear factorization scheme to the /cr-factorization scheme [fj . 
During the last years evidence has been collected that, for certain final states in the small-x 
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region, the use of this scheme provides a better description of HERA data than the traditional 
collinear factorization scheme: for the modeling of final states it is therefore important to 
have the off-shell photon-gluon fusion not only in LO but also in NLO. 

(a) (b) 

1 loop i tree tree i tree 





FIG. 1. Different contributions to the 7* impact factor at NLO. 

Finally, it is well-known that the leading order photon impact factor can be interpreted 
in terms of the photon wave function and a color dipole cross section ||@]. In recent years 
this color dipole picture has proven to be extremely useful for analyzing HERA data in the 
small-x region, in particular the low-Q 2 region of the F2 structure function and the diffractive 
final states 0-0]. The validity of the color dipole picture (or, maybe, its generalization to 
higher multipole terms) beyond leading order therefore presents an issue which has to be 
investigated. In NLO, the validity of the photon wave function picture is one of the results 
that should come out of the NLO photon impact factor calculation: in the interpretation of 
our results we will therefore pay particular attention to this aspect. 



II. KINEMATICS AND CALCULATIONAL TECHNIQUES 

The notations of the scattering process are depicted in Fig ^. The incoming momenta 
q = q' — xp and p are used to construct the Sudakov basis q',p (with q' = p 2 = 0) for the 
other momenta, and s = 2p ■ q' denotes the center of mass energy of the incoming particles. 
The virtuality of the photon is q 2 = —Q 2 , and as usual, x = Q 2 / s. The momentum transfer 
is denoted by r (with r 2 = t), the outgoing gluon carries momentum P. and the momenta of 
the quark and antiquark are k + r and q — k — £, respectively. We express these momenta 
in terms of their Sudakov variables, 

k = aq + [3p + k± , (1) 

£ = a e q' + (3 e p + t± , (2) 

r = a r q + (3 r p + r± , (3) 

where the Euclidean notation is used for the two dimensional transverse momenta, e.g. 
k 2 = —k\. We are interested in the limit of large energy s: s has to be much larger than 
Q 2 , k 2 , £ 2 , and r 2 , and the a variables of the ggg-system are of the order unity. The regions 
of integration are 



0(t/s) < a < 1 



\£\ 



< olo < 1 — a. 



(4) 



q — k — £ q 




FIG. 2. Kinematics of the process 7* + q — > qqg + q. 



For the produced gluon we require a large rapidity gap between the gluon and the lower 
quark. In other words, we include both the fragmentation region of the virtual photon 

aio < ct£ < 1 — a (5) 

(where aio denotes a small but energy independent cutoff parameter to be specified later), 
and the 'upper' part of the central region 

—= < a e < a eo . (6) 

There is a remaining part of the phase space which we do not consider: the 'lower' half of 
the central region: 

£ 2 III 

7 <at< ^ (7) 

and the fragmentation region of the lower quark. The central region belongs to the LO 
BFKL calculation, whereas the fragmentation region of the lower quark is part of the quark 
impact factor ||. The definition of the 'upper' and 'lower' parts of the central region is 
somewhat arbitrary. We have chosen the 'symmetric' separation: ai = -y= corresponds to 
the symmetric point M 2 = (q + r) 2 = Mj = (p + £ — r) 2 = £ 2 s. 

From the on-shell conditions of the outgoing quarks and gluons we obtain a r = t/s and 

£ 2 
Pts = — , (8) 

Ps = - } k + t? ~--Q\ (9) 

1 — a — ai Oi£ 

a 1 — a — ot(. Oil 

The Feynman diagrams are listed in Fig. |3|. The graphs in the second line (denoted by 
B\, . . . , B§) are obtained from the diagrams of the first line (Ax, . . . , A§) by interchanging 
the quark and antiquark, i.e. by transforming k + r^q — k — £. In terms of the Sudakov 
components this transformation reads: 



A 1 




FIG. 3. Feynman diagrams for the process 7* + q — > qqg + q. 



a — > 1 — a — at (11) 

k^-(k + r + £). (12) 

In correspondence with this transformation we introduce the shorthand notation: 

Oil = a , (13) 

on = {1-a-ai) , (14) 

a 2 = (1 — a) , (15) 

a 2 = (a + at). (16) 

where the bars stand for the transformation (|TTj), (|T2"D. The denominators, that will occur 
during the calculation, have the form: 



Di = (k + i + r-q) 2 

( 
D 2 = {k + r-qf 



- ai Q 2 -(k + £ + r) 2 -^(k + r) 2 -^£ 2 



Oil 



Q( 



= -a 2 Q 2 - — (k + r) 2 
D 3 = (k + £ + r) 2 

= -(k + £ + r) 2 + ^(k + r) 2 + ^£ 2 



(17) 



(18) 



Oil 



Q( 



a^ai 



-a i (k + r) + ai£f 



(19) 



D, = (k- qf 



k 2 + ^(k + £) 2 + ^£ 2 
ai a e 



D 5 = (£- rf 



1 (a e (k + £) + a 1 £) 2 (20) 



, Qt> - n ( fc + r) 2 _ n ( fc + £)» _ (/ _ r) 2 ( 2i) 



D B = fc 2 



12 r.2 "1 /i. , e\2 "1 ,2 



= _ ai Q2_ fe 2_- (fe + £) 2_- £ 2 (22) 

£> 7 =(A; + £) 2 

= -a 2 Q 2 - — (fc + £) 2 . (23) 

Under the transformation (p]) -([l2]) the denominators D\ to D-j obey the following relations: 

Di = D 6 , (24) 

^2 = D 7 , (25) 

D 3 = D 4 , (26) 

A, = ^s • (27) 

Rather than beginning with the scattering amplitude we immediately turn to the squared 
matrix elements, averaged (summed) over the helicities and color of incoming (outgoing) 
quarks. The polarization vectors of the photon are 

s L = ^(q' + xp), (28) 

e± = -^(0,l,±i,0). (29) 

In fact, since El = q/Q + 2xp/Q, we make use of gauge invariance and simply put El = 
2xp/Q. In this paper we restrict ourselves to the longitudinal polarization; the transverse 
case will be presented in a subsequent paper. 

Our calculations are performed in the Feynman gauge. As to the internal (t-channel) 
gluon numerators, we use the decomposition 

2 
0/u/ = - (Pfiql + ?»%) + V ( 30 ) 

In the high energy limit we retain only the first term, since the contributions of the other 
terms will be suppressed by powers of t/s. For the diagrams where the produced gluon 
couples to the lower fermion line, this may not be immediately obvious, since some of the light 
cone components of the t channel gluon are not small; we have checked the approximation 
against a calculation involving the full numerator, and we have obtained the same answer. 
In the diagrams A 5 and B 5 we have combined two Feynman graphs (Fig. |J) into one effective 
diagram. Since the produced gluon is far from the lower quark in rapidity, we can simplify 
the emission vertex: 



r — £, a 





Lc,\ 



P ^^ p — r 

FIG. 4. Emission of a gluon from the lower quark line. 



u(p — r, h!) 



'M - 



-t a t c 



7 A (^ + / 



rt a 



u(p, h) 



2p* 



[t c ,t a }S, 



hi,' 



(31) 



(p-£) 2 (p + £-ry at 

Here we have used that I has large components only in the q' direction; h and h! denote the 
helicities of the lower incoming and outgoing quarks, respectively, and the Lorentz index A 
is to be contracted with the polarization vector of the outgoing gluon with color index c. t a 
and t c are the color group generators in the fundamental representation, and a denotes the 
color index carried by the t-channel gluon. 

The result for the helicity-summed square of the matrix element (averaged over helicity 
and color of the incoming quark) has the following form: 

5 1 n 2 S ab 

\M\ 2 = e 2 e)g 4 5 ab 16 Q 2 s 2 ^ {AA VJ + AB l3 + BA t] + BB tj ) - A |^-, (32) 



where eef is the charge of quark flavour /, and g is the strong coupling constant. Following 
the labeling of Fig. [| we denote the product of graph Ai with graph A, by AAij etc. 
Furthermore, the phase space over the final state (including the ^-functions from the on- 
shell conditions) is given by: 

dada e d D ~ 2 kd D - 2 £d D - 2 r 



d$A 



1s(2<k] 



3D-4 



OL\OL\OL£ 



(33) 



D denotes the space time dimension. We calculate our results in dimensional regularization 
with transverse dimension D — 2 = 2 — 2e. 

The sum on the rhs of fl32"D , together with the phase space ( |3~3"1) and the flux factor l/(2s) 
represents the (helicity-summed) cross section of the process 7* + q — > (qqg) + q: 

d a = —\M\ 2 d§ 4 . (34) 

The NLO real corrections to the photon impact factor will be obtained from the sum on the 
rhs of ( |3"2"D by doing the integral over the qqg phase space and invariant mass. However, we 
will have to subtract the production of the gluon in the central region (multiregge kinemat- 
ics). This piece belongs to the LO BFKL calculation and must not be counted twice (see 
below). Integrating over the momentum of the outgoing gluon leads to infrared singularities 
which cancel against the infrared singularities of the virtual corrections || . Since the calcu- 
lation of the virtual corrections has been done using the dimensional regularization, we will 
use these same techniques also for the real corrections. 

In order to obtain the NLO jet cross section for the off-shell photon-gluon fusion process 
7* -\-g — > (2 jets, 3 jets), we have to combine the sum in (|32|) with the virtual NLO corrections 
to 7* + g — > qq, using a suitable jet definition. 
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III. RESULTS 

In this section we list our results for the helicity-summed matrix element. The traces have 
been calculated with the help of the Mathematica package FeynCalc flO ]. The results have 



then been simplified further by casting them into a factorizing form: in order to have the 
photon wave function interpretation we have tried to find a representation in which each 
term AAij can be written as a sum of products of two terms: one factor belongs to the 
incoming photon, the other one to the outgoing photon, and in between we have a piece 
which belongs to the cross section which describes the interaction of the ggg-system with 
the lower quark. After some algebra we arrive at the following results: 

2 
AA U = C F -J- 2 (1 - e) a x on A 2 (35) 

IJ 1 L> 2 

A4i2 = - — 2 ( (1 - e) ocx a 2 A x A 2 - a x a\ a x r 2 ) (36) 

iV c L)\ L) 2 J-^4 

A4i 3 = C F n2 ^ - ((l-e)a 2 1 A 1 A 3 + a\ (k + r) (k + £ + r)) (37) 

JJ 1 L> 2 L>3 

AA U = - N c ( - — a x (a 2 ax + 2 (1 - e) of) A x A 4 

H — ax a 2 a\ D 2 H ax a\ ax A x r — ax a\ ax r 2 ) (38) 

2 2a e J 

AA 15 = N c ax a\ a\ r 2 (39) 

at Dx D L 2 D 5 

*** = C F „_ n 2 2n2 (1 - e) ax a 2 2 A\ (40) 



, , ,, --,m,r" 2 " 2 

olxD\D\ 



AA 23 = -— ((1 - e) a 2 ax A 2 A 3 + a 2 «i {k + r -air) (k + £ + ax r)) 

l\ c L>x L> 2 L>3 L>a 



(41) 



AA24 = N c 2 ( — ax a 2 (a 2 ax + 2 (1 - e) a 2 ) A 2 A 4 

D?,DxD 5 \2a e ax v ' 

H — «! a 2 ai _D 2 H «i «2 <^i 4-2?" + «i «2 ^i r2 ) (42) 

2 2a f / 

-4-425 = - iV c — n2 n n ai a\ ax r 2 (43) 

a* D{ D A D 5 



•M* = C F - n2n2 (1 - e) a? A 2 (44) 



ai D{ D 



2 n2 

3 



1 / 1 _ 

A4 34 = iV c - — ai (ai ai - a £ - 2 (1 - e) a 2 ) A 3 A A 

DxL> 2 D 3 D 5 \2a e 

O A. r\ 

— —OLiQ-i a 2 D 2 + a x a 2 a ± D 7 + - — a x 0:2 a t A 3 r + ax a 2 a l r 

2 2ai 

+ ^a 2 axA 3 (k + eA (45) 



AA 



:-!.-> 



N r 



1 



AAaa = N. 



a e D 1 D 2 D 3 D 5 
1 / 2 



0L\ a 2 &l ®- 2 T 



D\D\ 



I — A 2 . «! ((1 — e) a? + a 2 «i) + «i «2 on ot\ D 2 + 2 a^ a, a i ^ 2 I 



AA45 = - N, 

AA 55 = 



c JyTJyi "i a l "i r2 



(46) 

(47) 

(48) 
(49) 



AB 12 
AB 13 
AB 14 



1 



1 



N c D 1 D 2 D 6 D 7 

-c F 2 



((1 — e) a\ ol\ A\ Bi + a.\ ol\ ik + a± r) (k + £ + a 2 r)) (50) 



1 - e) «i a 2 A 1 B 2 + «i a 2 (fe + r) (k + £ + r)) 



1 



£>i D 2 D 3 D 7 

( (1 — e) « 2 Ax £? 3 — a\ a 2 oi\ r 2 ) 



N c D x D 2 D A D 



N r 



n n n n ( T> — «i («i «i - «/ - (1 - e) 2 a 2 ) A x B 4 
D 1 D 2 D 5 D 7 \2a>e 



+ «i a\ a 2 D 2 «i a: 2 tt i -^7 H a i a 2 ol\ a 2 A\T — a.\ a 2 a.\ a 2 r 



2 a t 



AB 15 
AB 22 
AB 23 

AB 24 



-a 1 a 2 A 1 (k + r) 



N r 



-a\ a 2 a\ a 2 r 2 



(51) 
(52) 



(53) 
(54) 



a e D 1 D 2 D 5 D 7 

~W r> r> r, r, d 1 ~ e ) a<2 " 2 A<1 B 2 + Oi 2 a 2 (k + r - a x r) (k + £ + a x r)) (55) 
iV c L> 2 L> 3 L) A L> 7 

2 , 

C F „ n - (l-e)a(a 2 A 2 B 



ai D 2 D\ D, 



(56) 



4-^6 



-N, 



^2^4 



1 / 1 

n n ~ — a 2 («i ai - a/ - 2 (1 - e) a 2 ) A 2 B 4 



1 



+ «! a 2 «i «2 -D2 ~~ 77CK1 a 2 a i -^7 + 7; — a i a 2 «i «2 A 2 r" + ot\ a 2 ot\ a 2 r 



2a e 



- -a 2 a 2 A 2 (k + r) 



AB 25 
AB 33 
AB U 



N r 



1 



-«! a 2 «! a 2 r 



(57) 
(58) 



atD 2 D 4 D 5 D 7 

1 1 

-77- n n n n ((1 - e) ai ai A 3 B 3 + a x a : (fe + r - ol\ r) (k + £ + a x r)) (59) 

iV c L>i L> 3 L> A JJq 
1 



~N r 



1 



-a 



2 (ai a 2 + 2 (1 - e) a 2 ) A 3 B 4 



AB35 = iV c 



D ± D 3 D 5 D 7 \2 ai a e 

1 1 

-ai a 2 a 2 -D 7 H ai a 2 a 2 A 3 r + « x a 2 a 2 r 2 

2 2«; 

— 2 —2 2 

-«i a x a 2 7* 



D l D 3 D b D, 



(60) 
(61) 



AB 4A = -N c ( (2 at ai + a e - 2 (1 - e) a 2 ) A 4 B 4 

L) 2 JJ 5 L>7 \ 



I I 

+ — oti a^ ai a2 -D2 + 77^1 ol 2 a^ «i -D7 + 2 a 4 a 2 «i 02 r 2 — a 2 (k + £) (fe + r) 



(62) 



Ai3 45 = N c 2 ai a 2 ai a 2 r 2 (63) 

1^2 U h L> 7 

AB 55 = . (64) 

Here we have used the following abbreviations: 

A 1 =a 2 £ + a e (k + r) (65) 

A 2 =a 1 £ + a e (k + £) (66) 

A 3 =a l £-a e (k + r) (67) 

A 4 = a 2 (€ - r) + a* (fe + r) (68) 

Bi=a 2 £-«/(& + £) (69) 

B 2 =ai£-a^(fc + r) (70) 

B 3 =a 1 £ + a e (k + £) (71) 

B 4 = a 2 {£ - r) - a e {k + £) . (72) 

In the above definitions, the vectors Bj are obtained from the A, by the transformations 
(PI), (0)- Note also the additional relations 

A 2 = B 3 , A 3 = B 2 . (73) 

The remaining terms AAij (with i > j) and jBAjj (with i > j) follow from the symmetry 
properties AAij = AAji and A£>ij = BAj%\ in order to find BAij and ABji (with i < j), and 
all the BBij's we use the substitutions (Pf), ([12]) . 

We finally note that in the sum of the four (zj) = (44) matrix elements: AA 44 + A6 44 + 
BAu + BB 44 , the term proportional to D 2 in A4 44 , the terms proportional to D 2 and .D7 
in AB44, and the analogous terms in £>A 44 and £>£> 44 cancel against each other. 



An important feature of the results (35)— (|64) is the ultraviolet behavior of transverse 



momenta: individual terms, namely AAn and AA 22 , AA\ 2 and AA21, and AAu, AA 24 , 
AAn, AA 42 , and AA 44 go as l/£ 2 for £ 2 ^> A; 2 . When integrating over £, they lead to 
logarithmic divergencies. However, in the sum these divergencies cancel, and all transverse 
momentum integrals are ultraviolet finite. 

IV. THE CENTRAL REGION AND THE ENERGY SCALE 

As we have discussed in the beginning of section II, our results are valid in the fragmentation 
region of the photon and in the 'upper' part of the central region. Before we move on to 
any further evaluation of our results we have to address the question of how to remove the 



central region. Let us, therefore, first study the limit where the emitted gluon belongs to 
the central rapidity region: 



\£\ 



<ae<a i0 



P n \£\ 

sa £ o V s 



(74) 



In this limit we should find agreement with the leading logarithmic approximation. First 



of all, we note that the phase space integral in fl33|) contains a factor 1/af. this makes 



the at integral divergent, since, as we will see, some of the AAij have a finite but non 



vanishing limit for at — > 0. Let us see how this works in detail. In the limit of small at, the 
denominators D 1; . . . , D 7 ((fi7|)-(|23|)) are expressed in terms of the two denominators D 2 and 
Dj; they in turn are closely related to the denominators D(k + r) = a(l — a)Q 2 + (fe + r) 2 
and D(k + £) = a(l — a)Q 2 + (k + £) 2 , that are well known from the Born approximation: 



Di 

D 2 

D 7 



I — a 



D(k + r) 

a 

-£\ 

an 

1 -^£ 2 , 
ai 



it-r) 



a 



on 



1 — a 



D(k + £) 



(75) 



In the numerators on the rhs of (35)-(]64]) we expand in powers of at and keep only the 
leading term. The contributions from the diagrams containing A 3 or B 3 drop out completely, 
since they involve an extra fermion propagator that gives a suppression by a power of ag. 
Similarly the diagrams with Ai, A 2 (or B\, B 2 ) on both sides. Nonzero contributions come 
from diagrams involving A± or B^, but they cancel in the sum. We are finally left with the 
contributions coming from A5 and B$. As to the »4v4ij- terms, all non vanishing terms have 



the common factor , n vi ; 

rhs of (|3^) we find the simplified form 



all non vanishing ABij terms contain 



d 2 d 7 - 



For the sum on the 



I .Ml 



2s 2 2(2tt) 



8 



(lb 



2N r . 



-N r 5 



(lb 



a(l-a)\q L (k + £,a)\ 2 J- 



(£ - r)- 






(76) 



Here, 



\y L (k + £,a)\ 2 = 6 -^a 2 (l-a) 2 Q 2 



7T- 



D(k + r) D(k + £) 



(77) 
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is the square of the leading order (longitudinal) photon wave function. The matrix re- 
does not discriminate between the coupling of the t-channel gluon to the upper quark or 
antiquark lines, i.e. the difference between Ai and Bi is contained entirely inside the photon 
wave function. The matrix itself has the rather simple form 



/0 













_1 

2 

r! 



2 + 2 



r 2 



{t-rf 



The sum of all diagrams is then given by 



\ 







(t-rf 



/ 



(78) 



|Mp = ^ s ^ 



a(\ 



a) 



\V L (k + l,a) 



V 



AN c r 2 (£-r) 
I 2 



1 g 2 5 ab 

r 1 2N r 



(79) 



^l\ 2 , including the color factor 5 ab gives 



This is exactly what we expect: the factor a(l — a 

the impact factor for virtual photons in the LLA. The factors l/(r — £) A and 1/r 4 are the 

propagators of the reggeized gluons in the t-channel, and 4N c r 2 (r — £) 2 /£ 2 is the square 



of the Lipatov vertex ||11|| . Finally, 5 ab /(2N c ) is the LLA quark impact factor (averaged 
(summed) over initial (final) colors and helicities). Together with the an integral in the 
region (|74]) which gives a ln-s factor, we recover the LO BFKL result. This demonstrates 
that our results correctly reproduce the central region of the produced gluon. It may be 
interesting to note that (|79|) may also be written in the factorized form: 



\M\ 



__ 2s 2^f g A 5 a b a{l 



a 



x4N r 



_(r-ey 



£_ 



V L (k + £,a)\ 2 
r-£ 



.(r-iy 



£_ 



1 g 2 5 ah 

r 1 2N r 



(80) 



M; a product of two terms, one belonging 



In this way it has the same structure as (p5|)-( 
to the incoming photon, one to the outgoing. 

In order to avoid double counting we have to remove this central region from our results. 
Details of this subtraction procedure will be presented elsewhere, here we only sketch the 
general argument [fl2|JT3[1 . Consider, as an example, the contribution AA15 which, as an -^ 0, 
has a non vanishing limit. The integral over an diverges logarithmically, but the Ins-piece 
belongs to the LO BFKL calculation. We therefore decompose: 



1-0 



— AAi 5 (a e ) 



l-a 



21 



da 1 

a t 



[AA 15 (a e ) - AA 15 (0)6(a eo - a e )] + A4 15 (0) ln(V£ 



Oil o - 

\£\ ' 



51) 



where ai denotes an energy independent scale much smaller than unity which may depend 
upon the other momenta, e.g. £. As an example, we could define, as the boundary between 
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the fragmentation region of the photon and the central region, a cutoff on the invariant mass 
of the qqg system, M$, which we choose to be larger than Q 2 , I 2 , k 2 and r 2 . In this case 



Oil o ~ T7i- I n the first term on the rhs, we can safely put to zero the lower limit of integration 



M 







(the error is of order 1/a/s), the second term is part of the LO BFKL calculation, now with 
a definite scale inside the logarithm. This example illustrates the role of the energy scale 
within the result of the NLO calculation: a change in the scale results in a finite term. As 
an example, starting from (|8T| ) and adding the analogous contribution from the 'lower' part 
of the central region (where we have an analogous cutoff parameter fyo), we can write: 

i i r aeo \ , i i rfo°\ i t s N,i A r \\ r -£\ a toPeos , co x 

(v ^T^ (v ^T^ = Vlk-V ( 1 2 )- ( } 

From this we see that our NLO correction does not fix the energy scale, but, for a given 
choice of the scale, it provides a well-defined answer. Moreover, it then determines how this 
answer changes under a modification of the energy scale. 

This feature becomes more transparent if we generalize to infinite order in a s . Let us 
consider, instead of the process 7*g — > (qqg)q, the scattering of two virtual photons with 
virtualities Q\ and Q\, i.e. we replace the lower quark line by another photon impact factor. 
In the LO approximation, the scale of the energy factor in the total cross section is not fixed, 
i.e. any change of the scale exceeds the LO accuracy. One possible choice is the symmetric 
scale So = kk', k = \k\, k' = \k'\: 



CX^**v* ( s ) 



5+io ° duj fd 2 k fd 2 k' 



"> 7 ^ (2n) 2 J s _ zoc 2m J k 2 J k> 2 

x (-^Y %*(k,Qi)Gu.(k,k')® r (k',Q 2 ). (83) 

This choice seems to be most natural from the point of view of the Green's function. On the 
other hand, for a high energy scattering amplitude it is more natural to chose the 'physical' 
scale so = \JQ\Q\- In NLO, a change of the scale affects both the impact factors and 
the Green's function, i.e. the definition of these quantities becomes scale dependent. For 
example, starting from ( |8"B| ) where all elements are defined with respect to the scale (kk'), 
and changing to the 'physical' scale \/Q\Q\ we obtain: 

n5+lco diu fd 2 k fd 2 k' 
a, 



" w (2*r)* A-*. 2wiJ V J k" 

x iokJ^ ikM (£)"<wk,*)M*,«>> (f )" • m 

It seems natural to absorb the additional factors (^-) w and (^f) w into the impact factors of 
the photon. We introduce the Mellin transforms 

^W^aKf) 7 ^' (86) 
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where x{l) is the eigenvalue function of the BFKL kernel (for simplicity we ignore the scale 
dependence of the strong coupling a s ). Substituting these Mellin transforms in the cross 
section formula ( 8j ) an d evaluating the integrals over k and k' (which lead to (5-functions 



in the 7's), we obtain the following result for the cross section, 

f s 1 1 f s+ioo duj ( s Y 

^W=7^ra ^7 7TTT ( 87 ) 



^7 ^?V^(7-!-!M-7-!-i) 



x 



2vn VQ2/ W-X(7) 



The change of the energy scale therefore leads to shifts in the 7 variables [|T^,|T^] • 

An important consistency is the agreement with the renormalization group equation. If 
we consider the DIS limit Q\ ^> Q\, our result has to agree with the NLO result of the 
DGLAP scheme, restricted to the neighbourhood of the point u = (here n = uo + 1 is the 
moment index). It is easy to see that in the limit Q\ ^> Q\ our formula can be rewritten in 
the following way: 

This last equation provides a stringent test for the choice of the cutoff parameter of ag. 
Previously, in [fT2| , [3~3|1 the impact factors for on-shell quarks and gluons have been considered: 



here the cutoff parameter has been derived from the angular ordering condition, together 
with the requirement that the collinear singularities are reproduced correctly. 

V. THE ggg-COMPONENT OF THE PHOTON WAVE FUNCTION 

Let us next address the question how our results fit into the picture of the photon wave 
function and the color dipole cross section. We expect that, if this picture remains valid in 
NLO, our results should lead to a new Fock component of the photon, the qqg component. 
For a precise definition we first have to define what we mean, in NLO, by qq component. As 
an example, a qqg state of the photon where the gluon is either soft or collinear with the quark 
(or antiquark), should be counted as a part of the qq component. These contributions contain 
the infrared singularities which cancel against those of the virtual corrections. Consequently, 
the NLO definition of the qq component requires a 'resolution scale' which, in the qqg state 
of the photon, separates the qq component from the qqg component. 

As a first step, we consider our results (|35|)-(|64D, and ignore both the infrared singularities 
and the subtraction of the central region discussed in the previous section. We will show that 
our results can be expressed in terms of a new photon wave function, ip q q g (Q 2 ; pi, P2, a, a^), 
and the interaction cross section a q g g (a, ai, pi, p-z). To this end, we Fourier transform into 
transverse configuration space. We begin with the terms proportional to AiAj, which are 
present in all AAij (except for the elements with j — 5). As an illustration, consider the 
contribution AA\\ to (|3^) . Obviously we can write 

AA 11 = { 4 1 J? (25 mn C i ,a 1 a 1 ) ^. (89) 

JJ1-U2 -lA-^2 
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(where the sum over m and n is included). Defining the Fourier transform of the last factor: 

(2tt) 2 Dl D 2 ' [ j 

and using (|T7|) , flPSp, and (|35"D, we note that the dependence of the integrand jyjj- upon r 
is only through the combination k + r. We therefore change the integration variable k into 
fe' = k + r and obtain 



Q J J d kd " j p ^kjp 2 l {■"■!) n _j pir , (1) (r p. 



u a 2 a\ J J {2tt) 2 D t D 2 



where 



^ AQ - Pl , P2 , a , a ^N^JJ^e^^ (92) 

has no dependence upon r. Eq. ( f?2] ) therefore defines a new Fock component of the (longi- 
tudinal) photon wave function. N is the normalization constant which can be determined 
only in combination with the NLO qq Fock- component. With a similar result for the first 
term in fl89|), AA\\ takes the form: 



9JL AAn = f rf P' d ^ ^l^'e 1 ^^' n (Q 2 ; p[, p' 2 , a, a e )* (2S mn C F a 1 a 1 ) 

«2«1 J {^ ) 

x / ^^ e " plke ' ip2le ~ ipir 4lAQ 2 ^ Pi, p 2 , «, <*)■ (93) 



Integrating over the transverse momenta of the ggg system, we obtain ^-functions in the 
p- vectors which allow to do the p[, p 2 integrals: 

q2n2 J dadaedkdidr ^^ = J ^ ^ /"^ ^ ^W ^ V^ mn ^ ;n , (94) 
where 



(ii) i 



Pi, P 2 ) = ^ / rfW«' (25 m „^i) e— . (95) 



The final definition of a q g g . mn (a,ai, pi, p 2 ) requires an overall constant which cannot be 
determined at this stage. The incoming and the outgoing photons have the same transverse 
coordinates, i.e. the transverse coordinates of the qqg system remain 'frozen' during the 
interaction with the lower quark. We consider this as an important feature of the real 
corrections. 
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It is easy to see that a similar analysis applies to the terms proportional to AiAj in all 
other AAij. In summary, we transform the transverse momenta in the following way: 



_A ± 
D,D 



k' = k + r £' = £ (96) 



2 



Al k' = k + r £> = £--^- r (97) 



D 2 D 4 I - a 

D t D 3 ' 

A 4 
D 2 D 5 ' 



k' = k + r £' = £ (98) 

k' = k + r £' = £-r . (99) 



In this way we eliminate the dependence upon r from the Fourier transforms and obtain 
expressions which, after insertion into (|32j), can be interpreted as photon wave function 



components ip q q g - The r dependence is entirely carried by the interaction cross sections 
<J l qqg, similar as in the qq case. We carry out the transformation of the elements Bi in a 
similar fashion. The shifts in the k and £ integrations are just transformed according to 

©,©• 

As the next class of terms, consider the terms proportional to r 2 . Starting with AA\ 2 , 

we proceed as before, but the Fourier transformation is done only for D 1 D and D D , re- 
spectively: using the same shift of momenta as before, we again find the dependence upon r 
concentrated in phase factors, and the remaining integrals define new (scalar) contributions 
to the photon wave functions ip q J g and ip„J g which do not depend upon r. The factor r 2 goes 
into the cross section a qq - g ' . 

Next those contributions to AAij which contain D 2 or D 7 in the numerator: the first 
appearance is in ^L4i4, and a quick survey shows that these terms show up only if i = 4 or 
j = 4. Therefore, these pieces have to be part of the wave function component tpqqg'- e -g- 
jfn or j^ 7 D . Again, the dependence upon r can be factored out and leads to a phase factor 
which becomes part of the interaction cross section. 

Finally, we have contributions proportional to A^r which appear only for j = 4: following 
again the same steps as before, we take the vector A4 as a part of the wave function ipqq g , and 
we absorb r into the cross section o"^ 4 -*. The contributions proportional to scalar products 
such as (k + r)(k+£) are treated in a similar way: one of the two vectors belongs to the wave 
function (f//*))* ,the other into ^ (there is some freedom of how this division is done). In 
any case, the dependence upon r can entirely be absorbed into the interaction cross sections 

&qqg- 

As result, all terms appearing on the rhs of AAij, ABij etc. can be cast either into the 
form fl9"4]), or 

dadai / d 2 p 1 d 2 p 2 (V'SU*) <r&lrrfi>qq g ( 10 °) 



or 



f dadai fs Pl d^(i,%)^%. (101) 
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In other words, we have both scalar and vector-like contributions to the new photon wave 
function ip q q g . 

At this stage, the matrix of interaction cross sections a^' (here the indices i and j have 
to be extended in order to include AA, AB, BA, and BB, i.e. they run from 1 to 10) is not 
diagonal: so far our choice of the basis of the wave function components had been dictated 
by the Feynman diagrams from which we have started our calculations. Since the matrix is 
symmetric, we can always find a more physical basis in which the interactions are diagonal. 
In this basis, our cross section can be written as a sum of squares, each of which supports 
the wave function interpretation. In particular, the transverse distances in the incoming 
and outgoing wave functions are the same: they remain 'frozen' during the interaction with 
the quark target. Note that the definition of eigenstates of the qqg-system is determined 
by the interaction cross section: it is the interaction with the quark which filters out the 
eigenstates of the qqg- system. 

As we have indicated before, this analysis is not yet complete. To proceed further, we 
have to combine the results of this paper with the virtual corrections obtained in []]]: a 
suitable definition of the NLO qq state requires a subtraction of a part of the phase space 
of the qqg final state. Together with the subtraction of the central region, this may very 
well lead to a 'soft melting' of the 'frozen' transverse distances. This would not be totally 
unexpected: a soft gluon occupies a large region in transverse space, and it may take the 
form of a non-local cloud around the quark or antiquark. 



VI. CONCLUSIONS 

In this paper we have computed, for a longitudinally polarized virtual photon, the 
helicity-summed squared matrix element of the process 7* + q — > (qqg) + q. The calcu- 
lation is valid in the Regge limit with a large rapidity gap between the qqg-system. and the 
other outgoing quark. The results are of interest for several different lines of research. In 
each case, further calculations are needed which will be based upon the results of this paper. 

First of all, our results present the real corrections to the NLO corrections of the photon 
impact factor. In order to complete the NLO calculations we have to integrate over the 
final momenta of the ggg-system and to combine with the virtual corrections obtained in 
|jlj] . In the sum of the real and virtual corrections the infrared singularities have to cancel; 
furthermore, we have to subtract the contribution of the central region. 

Next, the results of this paper contain the three-parton final state of the photon-gluon 
fusion process 7* + g — > qqg where the incoming gluon is off-shell. The complete NLO 
calculation of this process also requires, in addition to the tree-level calculation of the process 
7* + g — » qqg, to consider the sum with the virtual corrections to 7* + g — ► qq (which is 
contained in Q). From the qqg final states we have to subtract the central region. With 
a suitable jet definition, one obtains the two-jet and the three-jet cross sections in the k?- 
factorization scheme. 

Finally, our results are part of the NLO corrections to the photon wave function picture. 
Ignoring infrared singularities and the problem of the central region, we have calculated the 
new qqg Fock component of the photon wave function. However, for a proper treatment of 
the infrared part we need to consider also the virtual corrections in |ff|: we have to define a 
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NLO qq Fock state which includes a part of the soft gluons in our qqg final state. 

Note added: Recently, a paper by V. S. Fadin, D. Y. Ivanov and M. I. Kotsky has 
appeared which also addresses the NLO calculation of the virtual photon impact factor ( [hep 



|ph/0106099|) . It summarizes the status of the virtual corrections, repeating the content of 
hep-ph/0007119| , and it presents the matrix elements for the process 7* + q — > (qqg) + q. 
In our paper we have calculated the helicity sum of the squared matrix elements for this 
process. 
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